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Abstract 

We discuss the status of the pion distribution amplitude (DA) from analyzing 
the CLEO experimental data in the context of QCD sum-rule techniques and 
QCD perturbation theory at the NLO accuracy. The constraints extracted this 
way for the Gegenbauer coefficients a% and 04 exclude $cz at the 4er level, 
while <£ aS y is outside the 3a error ellipse. These data provide strong support 
for the type of endpoint-suppressed, double-humped pion DA we derived via 
QCD sum rules with nonlocal condensates and favor a value of the vacuum 
quark virtuality X 2 ~ 0.4 GeV 2 . This pion DA is in agreement with the E791 
data, though these experimental results should be viewed carefully and further 
confirmation is necessary for a more accurate judging of pion DAs from them. 



1. INTRODUCTION 

It is widely believed today that the nontrivial QCD vacuum plays an important role in understanding the 
analytic properties of hadron distribution amplitudes (DA) in terms of their quark and gluon degrees of 
freedom [ 1 2|. In fact, one can use QCD sum rules with nonlocal condensates [3| to connect dynamic 
properties of (light) mesons, like form factors and DAs, directly with the QCD vacuum. The classical 
example is the pion DA, which describes how the pion's longitudinal momentum is shared between 
its quark and antiquark constituents when probed at large momentum transfer Q 2 . First, a detailed 
knowledge of the pion DA is necessary in order to make precise calculations of "hard-scattering" form 
factors, like F% m (Q 2 ) and F 7r7 (Q 2 ) and compare the results with the experimental data. Second, the 
parton structure is interesting in its own right, providing insight into the nonperturbative hadron structure 
at large distances. There is a long history of determining the pion DA stalling from the asymptotic limit 
of perturbative QCD El El to QCD sum rules EE] ED HE] ED, to instanton-based models fTOl fTTI fT2l . 
and to lattice computations [13|. Now, with the help of the recent CLEO data [14|, one gets a handle 
on the pion DA from the experimental side. Indeed, these data can be analyzed fl31 IT6l using the 
framework of light cone sum rules [17| to extract constraints on the Gegenbauer coefficients 02 and 
of the conformal expansion. Supplementary constraints on the shape of the pion DA are supplied by 
diffractive di-jets production events [18|, but the uncertainties of these experimental results and their 
controversial theoretical interpretations Q9]|20]|2] are still prohibiting definite conclusions. 



2. NONLOCAL CONDENSATES: THE MODUS OPERANDI 
2.1 Modelling the nonperturbative QCD vacuum 

The nonlocal quark condensate represents a partial resummation of the OPE to all orders in terms of the 
vacuum expectation value of the nonlocal operator 
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where a, p are spinor indices and the integral in the Fock-Sch winger string E(0, z) is taken along a 
straight-line path. Note that Fsy(z 2 ) are analytic functions around the origin and that their derivatives 



at zero are related to condensates of corresponding dimension. Recall that the condensates of lowest 
dimensions 
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form the basis of the standard QCD sum rules [ 22 ] and have been estimated, while higher-dimensional 
ones are yet unknown. In the chiral limit one has 
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and the parameter \ 2 /2 fixes the width of Fs(z 2 ) around the origin. 

For not too large Euclidean distances z 2 = — z\ > 0, the nonlocality behavior of the (quark) 
condensate can be implemented by the Gaussian ansatz [3| F^(z 2 ) = exp (— \ 2 z 2 /8) , in which X q has 
the meaning of an inverse vacuum quark correlation length. This corresponds to the specific form of the 
virtuality distributions / s = 5 (a - X 2 /2) , / v ~ a s Q 3 5' (a-X 2 /2) EHU with 
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This kind of virtuality distributions fix only one main property of the nonperturbative vacuum — quarks 
can flow through the vacuum with a nonzero momentum k, and the average virtuality of such vacuum 
quarks is just (A; 2 ) = A 2 /2 (for a determination of \ 2 „/2 from lattice data, see l23l ). 



2.2 Nonlocal QCD sum rules and the pion distribution amplitude 

The pion DA of twist 2, if n (x, n 2 ) (x being here the longitudinal momentum fraction), defined by 
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can be related to the nonlocal condensates by means of the following sum rule 
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where the index i runs over all scalar, vector, and tensor condensates with dim=6 |9] |S) ; M 2 is the 
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Fig. 1: (a) First nonzero moments {£ (up to N = 10, see text) and (x 1 )* = (l/3)(a; 1 ) 7r — 1 of i^bms (dark bars) with 
the upper and lower error-bars (grey bars), determined with nonlocal QCD sum rules |9|. (b) Histo gram of the first nonzero 
Gegenbauer coefficients a„ of the BMS pion DA and the envelopes of the "bunch" like in (a). 



Borel parameter, the duality interval in the axial channel. This sum rule allows us to determine 




Fig. 2: (a) Bunch of pion DAs for Aq = 0.4 GeV 2 determined with nonlocal QCD sum rules at normalization scale /i « 
1.0 GeV. (b) The BMS model |9 (solid line) in comparison with the asymptotic DA (dotted line) and the CZ model (dashed 
line) Q. 



the first ten moments {£, N ) n = J fn(x)(2x — l) N dx of the pion DA and independently the inverse 
moment (x -1 )^ = ip 7T (x)x~ 1 dx quite accurately (see in ll24l for an illustration). The corresponding 
Gegenbauer coefficients can be determined (see [9|) from this set within some error range (Fig.^>) that 
translates into a "bunch" of pion DAs shown for X 2 = 0.4 GeV 2 in Fig. El 

Their striking feature is that their endpoints are suppressed relative to both the CZ model and the 
asymptotic solution. Crudely speaking, this shape structure is the net result of the interplay between the 
perturbative contribution and the non-perturbative term A<3?s(x; M 2 ) (related to the scalar condensate) 
that dominates the RHS of the SR in Eq. ©. The fact that the function A<J>s(£; M 2 ) is not singular in 
x and has a dip at the central point of the interval [0, 1] is also reflected in the shapes of these DAs. We 
emphasize that a suppression of the endpoint region as strong as possible (for a dedicated discussion we 
refer to 12510 is important in order to improve the self-consistency of perturbative QCD in convoluting 
the pion DA with the specific hard-scattering amplitude for a particular exclusive process. In order to 
stress this point, we show in Fig.^x™ 1 )^, calculated as JJ + °'° 2 ip(x)x~ 1 dx and normalized to 100% 
(y-axis). The main message from this figure is that (x^ 1 )^ receives in the endpoint region, say between 




Fig. 3: Percentage distribution (see text) of the first inverse moment in x of the BMS model 1 9 1 in comparison with rival models. 



< x < 0.1, a contribution of only 17% in the case of the BMS model, whereas it reaches as much as 
40% for the CZ model and still 19% for the asymptotic solution. 



3. ANALYSIS OF THE CLEO DATA 

It was shown by Khodjamirian [17] that the light-cone QCD sum-rule (LCSR) method provides the 
possibility to treat the problem of the photon long-distance interaction (i.e., when a photon goes on 
mass shell) in the 7*(Q 2 )7(<? 2 ) — ► tt form factor by performing all calculations for sufficiently large 
q 2 , using quark-hadron duality in the vector channel, and then analytically continuing the results to 
the limit q 2 = 0. Schmedding and Yakovlev (SY) fTT) applied these LCSRs to the NLO of QCD 



perturbation theory. More recently, we have [16] taken up this sort of data processing (i) accounting 
for a correct Efremov-Radyushkin-Brodsky-Lepage (ERBL) fJJ |5j evolution of the pion DA to every 
measured momentum scale, (ii) estimating more precisely the contribution of the (next) twist-4 term, 
and (iii) improving the error estimates in determining the la- and 2<j-error contours in the (02,04) 
plane. Moreover, our error analysis takes into account the variation of the twist-4 contribution and treats 
the threshold effects in the running of a s {Q 2 ) more accurately. 

Our procedure is based upon LCSRs for the transition form factor F 7 * 77r (Q 2 , q 2 ~ 0) fT7llT5l : 
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following from a dispersion relation with M 2 m 0.7 GeV 2 , where m p is the p-meson mass and sq = 
1.5 GeV 2 denotes the effective threshold in the p-meson channel. The spectral density p(Q 2 , s; fi 2 ) 

1 11 1 Fq C q ^(Q 2 , q 2 = —s; 11 2 ) is calculated by virtue of the factorization theorem for the form factor at 
Euclidean photon virtualities q\ = —Q 2 < 0, q\ = —q 2 < 131 151 1261 . and the factorization scale u 2 is 
fixed by SY at u 2 = /x| Y = 5.76 GeV 2 . Moreover, F^^ n (Q 2 , q 2 ; p?) contains a twist-4 contribution, 
which is proportional to the coupling 5 2 (u 2 ), defined by fTTT [TTt (i:(p)\g s dG ail 'j o 'u\0) = ib 2 ^p p , 
where G aiX = (l/2)e appa G^ and G prj = G a pa X a /2. 

This contribution for the asymptotic twist-4 DAs of the pion as well as explicit expressions for 
the spectral density p(Q 2 , s; fi 2 ) in LO have been obtained in (T7i to which we refer for details. The 
spectral density of the twist-2 part in NLO has been calculated in [15| — see Eqs. (18) and (19) there. 
All needed expressions for the evaluation of Eq. © are collected in the Appendix E of [16|, cf. Eqs. 
(E.1)-(E.3). We set u? = Q 2 in F^^ 77 (Q 2 , q 2 ; p 2 ) and use the complete 2-loop expression for the form 
factor, absorbing the logarithms into the coupling constant and the pion DA evolution at the NLO level 

fl~6l so that ots^u 2 ) a s (Q 2 ) and (p^{x\ p 2 ) Pn( x > Q 2 ) = U(p? — > Q 2 )p n (x; fj?) (RG denotes 
the renormalization group). Then, we use the spectral density p(Q 2 , s,Q 2 ), derived in [ 15 ] at u? = /i| Y , 
in Eq. Q to obtain F 7 77r ((5 2 ) and fit the CLEO data over the probed momentum range, denoted by 
{Qexp}- m our recent analysis fl~6l the evolution ip 7T (x;Q 2 ) = U (p| Y — * Q 2 ) ( P-k(x; /x| y ) was performed 
for every individual point Q 2 xp , with the aim to return to the normalization scale /i| Y and to extract the 
DA parameters (02, 04) at this reference scale for the sake of comparison with the previous SY results 
fT5l . In effect, for every measurement, {Q 2 xp , F J 77r (Q 2 xp)}> its own factorization and renormalization 
scheme was used so that the NLO radiative corrections were taken into account in a complete way. 

The accuracy of this procedure is still limited because of the uncertainties entailed by the twist-4 
scale parameter 1151 . k ■ 5 2 , with the factor k expressing the deviation of the twist-4 DAs from their 
asymptotic shapes (another source of uncertainty originates from the unknown NNLO a s -corrections, 
see fToll ). Based on our experience with the twist-2 case, we set k = 1 ± 0.1. As a result, the final 
(rather conservative) accuracy estimate for the twist-4 scale parameter can be expressed in terms of 
k ■ S 2 (l) = 0.19 ± 0.04 GeV 2 fToll . To produce the complete 2a- and lcr-contours, corresponding to 
these uncertainties, we need to unite a number of regions, resulting from the processing of the CLEO 
data at different values of the scale parameter k ■ 5 2 within this admissible range [ 1 6 1 . The obtained results 
for the asymptotic DA (♦), the BMS model (K) fj), the CZ DA (■), the SY best-fit point (•) [T5J, a 
recent transverse lattice result (T) (fourth reference in fT3l ). and two instanton-based models, viz., (*) 
[10 1 and (♦) (using in this latter case m q = 325 MeV, n = 2, and A = 1 GeV) [11], are displayed 
in Fig. Ua) varying the twist-4 scale parameter k ■ 5 2 in the interval [0.15 < k ■ 5 2 < 0.23] GeV 2 . 
The important points to observe from this figure are these: (i) the nonlocal QCD sum-rule constraints, 
encoded in the slanted shaded rectangle, are in rather good overall agreement with the CLEO data at the 
ler-level; (ii) the CZ model and the asymptotic DA are ruled out at least at the 3<r-level. These findings 
are not significantly changed, even if one allows an extreme twist-4 uncertainty of 30%, or if excluding 
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Fig. 4: (a) Analysis of the CLEO data on F 1T1 » -y(Q 2 ) in terms of error regions around the best-fit point (+) (broken line: la; 
solid line: 2a; dashed-dotted line: 3a) in the (a,2,cn) plane contrasted with various theoretical models explained in the text. The 
slanted shaded rectangle represents the constraints on (a2, 04) plane posed by the nonlocal QCD sum rules |9| for the value 



A 



0.4 GeV . All constraints are evaluated at /i SY = 5.76 GeV after NLO ERBL evolution, (b) Comparison of tp ahy (solid 



line), ifi cz (dashed line), and the BMS "bunch" of pion DAs (strip, 1 16|) with the E791 data 1 18|. The corresponding \ 2 values 
are: 12.56— asy; 14.15— CZ; 10.96— BMS. 



the low-momentum-transfer data tail — say, up to Q^ xp = 3 GeV 2 1161 . In the first case, the asymptotic 
DA is outside the 3<r-error ellipse, whereas in the second case it remains outside the 2a region with the 
instanton-inspired models just at the 2<7-ellipse boundary and the CZ model always far outside. 



4. E791 DATA: CONSTRAINTS FROM DIFFRACTIVE DI-JETS PRODUCTION 

An independent source of experimental data to constraint the shape of the pion DA is provided by the 
E791 Fermilab experiment [ 18 1. Unfortunately, these data are affected by inherent uncertainties and their 
theoretical explanation by different groups fl"9"ll20ll271 is still controversial so that they cannot be used to 
exclude some optional model. For our exposition here the important point is to show that our predictions 
for this process are not conflicting the E791 data using for all considered models the same calculational 
framework, notably the convolution approach of ED - The results of the calculation are displayed in Fig. 
0Jb) making evident that the E791 data are relatively in good agreement with our prediction — especially, 
in the middle x region, where our DAs "bunch" has the largest uncertainties (see Fig.|2k). Note, however, 
that all theoretical predictions shown in this figure are not corrected for the detector acceptance. For a 
more precise comparison, this distortion must be taken into account. 



5. CONCLUSION 

Both analyzed experimental data sets (CLEO [ 14 1 and Fermilab E791 [18)) converge to the conclusion 
that the pion DA is not everywhere a convex function, like the asymptotic one, but has instead two 
maxima with the end points (x = 0, 1) strongly suppressed — in contrast to the CZ DA. These two key 
dynamical features of the DA are both controlled by the QCD vacuum inverse correlation length X q , 
whose value suggested by the CLEO data analysis is A 2 ~ 0.4 GeV 2 in good compliance with the QCD 
sum-rule estimates and lattice computations. 
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